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We present a numerical study of the electromagnetic re- 
sponse of the metamaterial elements that are used to con- 
struct materials with negative refractive index. For an array 
of split ring resonators (SRR) we find that the resonant be- 
havior of the effective magnetic permeability is accompanied 
by an anti-resonant behavior of the effective permittivity. In 
addition, the imaginary parts of the effective permittivity and 
permeability are opposite in sign. We also observe an iden- 
tical resonant versus anti-resonant frequency dependence of 
the effective materials parameters for a periodic array of thin 
metallic wires with cuts placed periodically along the length 
of the wire, with roles of the permittivity and permeability re- 
versed from the SRR case. We show in a simple manner that 
the finite unit cell size is responsible for the anti-resonant be- 
havior. 

PACS numbers: 41.20.Jb, 42.70.Qs, 73.20.Mf 

The recent development of metamaterials with negative re- 
fractive index - or double-negative (DNG) metamaterials 1 has 
confirmed that structures can be fabricated that can be in- 
terpreted as having both a negative effective permittivity e 
and a negative effective permeability fi simultaneously. Since 
the original microwave experiment of Smith et al., 2 various 
new samples were prepared, 3,4 all of which have been shown 
to exhibit a pass band in which the permittivity and per- 
meability are both negative. These materials have been used 
to demonstrate negative refraction of electromagnetic waves, 5 
a phenomenon predicted by Veselago. 6 Subsequent experi- 
ments 7 have reaffirmed the property of negative refraction, 
giving strong support to the interpretation that these meta- 
materials can be correctly described by negative permittivity 
and negative permeability. 8,9 

There is also an increasing amount of numerical work 10-1 2 
in which the transmission and reflection of electromagnetic 
wave is calculated for a finite length of metamaterial. For a 
finite slab of continuous material, the complex transmission 
and reflection coefficients are directly related to the refractive 
index n and impedance z associated with the slab, which can 
in turn be expressed in terms of permittivity e and perme- 
ability [i. A retrieval procedure can then be applied to find 
material parameters for a finite length of metamaterial, with 
the assumption that the material can be treated as continu- 
ous. A retrieval process was applied in Ref., 13 and confirmed 
that a medium composed of split ring resonators (SRRs) and 
wires could indeed be characterized by effective e and fi whose 
real parts were both negative over a finite frequency band, as 
was the real part of the refractive index n. 

The retrieval process, however, uncovers some unexpected 
effects. For the SRR medium, for instance, the real part of the 



effective permittivity e exhibits an anti-resonant frequency 
dependence in the same frequency region where the perme- 
ability undergoes its resonance. This anti-resonance can be 
seen in the composite SRR+wire negative index medium as 
well. The anti-resonance in the real part of the permittivity 
is also accompanied by anti-resonant behavior in the imagi- 
nary part of the permittivity e", which exhibits an absorption 
peak opposite in sign to that of the imaginary part of the per- 
meability. Assuming that waves have a time dependence of 
exp(— not), one would expect the imaginary parts of both the 
permittivity and permeability to have positive values at all 
frequencies, since the material is passive. 14, 15 

The aim of present paper is to show that the anti-resonant 
behavior of the material parameter is an intrinsic property of 
a metamaterial, a consequence of the finite spatial periodic- 
ity. To illustrate this point, we present the retrieved material 
parameters for two types of metamaterial media that have 
been used to form negative refractive index composites: the 
first medium comprises an array of SRRs that exhibit a reso- 
nant permeability and anti-resonant permittivity. The second 
medium comprises an array of cut metallic wires, which ex- 
hibits resonant permittivity and anti-resonant permeability. 

We used the transfer matrix method to simulate numeri- 
cally the transmission of the electro magnetic waves through 
the metamaterials. Transmission data for an array of SRR 
were published elsewhere 10, 12 and will not be repeated here. 
To the best of our knowledge, the analysis of the array of cut 
wires has not been presented in the literature yet. Fig. 1 
shows the frequency dependence of the transmission and ab- 
sorption for an array of cut wires. As expected, the system 
exhibits a band gap for frequencies fo < / < fp. Here, /o is 
the resonance frequency and f p is the plasma frequency. The 
transmission data are similar to that for an array of the SRR. 
However, contrary to an array of SRR, the system exhibits 
the resonant behavior of the effective permittivity at f = f . 

From the transmission and reflection data we calculate the 
effective permittivity and permeability. Details of the method 
were published elsewhere. 13 The method is based on the as- 
sumption that the system is homogeneous. Textbook formu- 
las for the transmission and reflection of the slab of width 
d are then inverted to obtain effective impedance z and ef- 
fective refractive index n. Permittivity and permeability are 
obtained from relations 

n — y/ejl and z — \J fi/e. (1) 

Typical frequency dependence of the effective parameters of 
an array of SRR and array of cut wires are shown in fig. 2 and 
3, respectively. Both structures exhibit qualitatively the same 



behavior: there is a resonant frequency interval, in which one 
effective parameter is negative (fjf for SRR, e' for cut wires). 
The resonant behavior of this parameter at the left border of 
the band gap is clearly visible. The frequency dependence of 
the second parameter is anti-resonant. Its real part decreases 
to zero at / = /o and imaginary part is negative for f > fo- 
Notice the qualitative similarity of both systems: data are 
qualitatively the same, and differ only in the exchange of e 
and /i (which is equivalent to the transformation z — > 1/z). 

The anti-resonant behavior of the effective parameter has 
its origin in the finite lattice period a associated with the 
metamaterial structure. One manifestation of the lattice pe- 
riodicity is that there is a maximal wave number, given by 
fcmax = ir/a. 16 If we assume that the metamaterial can be 
treated as a continuous medium with an index of refraction n, 
then the definition of n = cn s htk/u> shows that n is necessarily 
bounded. The generic ui(k) dispersion diagram of a resonant 
periodic structure results from the coupling of a dispersionless 
resonance curve at the resonant frequency uio = 27r/o and the 
light line / = cu s ht/k (see Ref. 2 for details). The result is 
a lower branch that extends from zero frequency to the reso- 
nance frequency fo , followed by a band gap that extends from 
fo to f p , followed by an upper branch that extends upwards 
in the frequency from f p . 

Because the resonant frequency of a typical metamaterial 
element implies a free space wavelength much longer than the 
unit cell size, an effective medium approach has been applied 
that results in a characterization of metamaterial in terms of 
bulk e and fi. The square root of the product e/i is the re- 
fractive index, which must be consistent with that determined 
from the dispersion diagram. To understand the generic prop- 
erties associated with a resonant system in a periodic struc- 
ture, we start by considering the dispersion characteristics of 
the system near the resonant frequency, where we have 
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where a is a real number. Solving Eq. 2 for k and using n — 
cughtfc/w, we find an approximate expression for the refractive 
index: 



n (f) ~ 7TT- ( - ~ oV^o - w ) 



(3) 



The maximum value of the refractive index at the resonance 
frequency is thus 
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determined only by the resonant frequency of the element and 
the periodicity. 

For the composite structures analyzed in this paper, Eq. 4 
gives n max = 4.24 for an array of SRRs and n max = 3 for the 
cut wires array. When comparing n max with data presented 
in figs. 2 and 3, we see that n' indeed does not exceed these 
limits. 

Assume that the effective parameter x (x represents the 
effective permittivity for the cut wires system and the effective 
permeability for the SRR system) exhibits a resonant form 
corresponding to 



Ff 2 



(5) 



where 7 is a damping factor, and F is a filling factor (fraction 
of volume of the metallic components). At resonance, the 
imaginary part x"(u>) becomes 



x"(u>) = (1 - 
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(6) 



(ojp = u>q/(1 — F)) and is greater than zero, as expected. 

If we require that the index n calculated form the disper- 
sion curves be consistent with the bulk permittivity and per- 
meability parameters, than we must have n(uj) = W x(u))y(u>). 
Near the resonance frequency this implies that the second ef- 
fective parameters y(ui) behaves as 



2/H = 



n 2 (cu) 
x(u) 



(7) 



Comparing the expression for y(uj) with that for x(u>) we see 
that the poles and zeros of x and y are reversed, as long as 
n is bounded with the form given by Eq. 4. Moreover, for 
small n"(uj), we immediately see that the product of imagi- 
nary parts x"(u>)y"(u)) < is negative. To be more specific, 
with the help of Eq. 1 we obtain that 
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One sees that the sign of e"/x" is fully determined by the r.h.s. 
of Eq. 8. We are not aware about any physical requirement 
which prevents the r.h.s. of Eq. 8 to be negative. In fact, 
data presented in Figs. 2 and 3 show clearly that 

\n"z'\ < \n'z"\. (9) 

in the left part of the resonance gap, since \z'\ < \z"\ and 
n" < \n'\. Thus, in agreement with Eq. 7 we conclude 
that the opposite sign of e" and fj," is a consequence of small 
transmission losses in the structure. The same conclusion was 
derived in Ref. for the DNG metamaterials. 

In conclusion, we presented the numerical analysis of the 
effective parameters of two metamaterials: an array of SRR 
and an array of thin metallic cut wires. We show that the 
effective parameters of these systems exhibit resonant and 
anti-resonant behavior similar to that found recently in the 
double negative metamaterials. We suggest that the anti- 
resonant behavior is caused by the requirement that the re- 
fractive index must be bounded in the structures which posses 
finite spatial periodicity. As the spatial periodicity is an un- 
avoidable property of the metamaterials, we conclude that 
observed seemingly unphysical behavior of effective material 
parameters is an intrinsic property of composites, which can 
not be avoided for instance by decreasing the size of the unit 
cell. 

The electromagnetic response of metamaterials is usu- 
ally embodied in a description involving bulk, continuous, 
frequency dependent permittivity and permeability tensors. 
This description, however, is only approximate, as spatial 
dispersion is always present to some degree in metamateri- 
als. Thus, the applicability of our retrieval procedure, which 



uses the formulas for transmission and reflection of the ho- 
mogeneous slab, might not be applicable in the neighbor of 
the resonance frequency, where the wave length of the elec- 
tromagnetic wave inside the composite is already comparable 
with the spatial period. 

Besides the analysis of the effective parameters, our results 
are interesting also for further development of new double 
negative metamaterials, in which the lattice of cut wires might 
find new interesting applications. 12 ' 18 
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FIG. 1. Transmission of the EM wave through the periodic 
lattice of thin metallic wires. Wires are parallel with the y 
axis, and the EM wave, polarized with E || y, propagates in 
the z direction. The system is infinite in the x and z direction, 
and 60 rows of wires are considered along the propagation 
direction. The structure is characterized by four length pa- 
rameters: the wire thickness w, gap A, the wire length L and 
the lattice period (mutual distance of wires) a. In the present 
simulations, the lattice constant is a — 3.66 mm, the thickness 
of the wire is w = 0.33 mm, the length of the wire is L = 7 
mm and the cut of the wire is A = 0.33 mm. Data show the 
drop of the transmission (solid line) at frequency fo w 13.35 
GHz. Dashed line is absorption. We found that the resonant 
frequency fo is almost independent on the lattice constant a 
(1.66 < a < 5 mm). However, it depends strongly on the gap 
A. In inset shows how the frequency fo depends on the ratio 
A/L. Solid line is the fit f = ^lg| [a ln(L/A) + aj] -172 
with ao = 0.48 and ai = 0.17. Plasma frequency f p depends 
on the lattice period a. For a = 3.66 mm we found numeri- 
cally f p « 24.5 GHz which agrees with predictions of Sarychev 
and Shalaev, 17 indicating that the value of the plasma fre- 
quency is not influenced by the wire cut. 
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FIG. 2. Effective parameters for an array of split ring res- 
onators (solid lines: real part, dashed lines: imaginary part). 
The resonant behavior of the effective permeability fx at fre- 
quency fo = 9.66 GHz is clearly visible. Shaded area shows 
the resonance frequency interval in which [i! is negative. Note 
the anti-resonant behavior of e. Note also that e" is nega- 
tive. The sharp discontinuity in e" is due to the extremely 
fast decrease of n' to zero in the right part of the resonance 
interval. The size of SRR is 3 mm, the size of unit cell is 
L x x L y x L z = 3.33 x 3.66 x 3.66 mm. The SRR lies in the 
x — plane, EM wave propagates along the z direction and 
is E || y polarized. 
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FIG. 3. Effective parameters of the periodic lattice of cut 
wires. Resonant behavior of the effective permittivity e as 
well as anti-resonant behavior of effective permeability [i are 
clearly visible. The real part of permeability p! is zero at fo, 
and the imaginary part of the permeability /i" is negative for 
/>/o- 



